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Abstract 

We investigate the time T a quantum computer requires to factorize a given 
number dependent on the number of bits L required to represent this number. We 
stress the fact that in most cases one has to take into account that the execution time 
of a single quantum gate is related to the decoherence time of the qubits that are 
involved in the computation. Although exhibited here only for special systems, this 
inter-dependence of decoherence and computation time seems to be a restriction 
in many current models for quantum computers and leads to the result that the 
computation time T scales much stronger with L than previously expected. 
PACS: 42.50.Lc 

I. Introduction 

Since Shor's discovery [|^, Q of an algorithm that allows thefactorization of a large number 
by a quantum computer in polynomial time instead of an exponential time as in classical 
computing, interest in the practical realization of a quantum computer has been much 
enhanced. Recent advances in the preparation and manipulation of single ions as well as 
the engineering of pre-selected cavity light fields have made quantum optics that field of 
physics which promises the first experimental realization of a quantum computer. Several 
proposals for possible experimental implementations have been made relying on nuclear 
spins, quantum dots IQ, cavity QED and on ions in linear traps p. 

One can estimate the time T needed for a single run of Shor's algorithm to be equal 
to the time t^i required to execute an elementary logical operation multiplied by the 
required number of elementary operations, which is of the form eL^ + 0(L^) 0. It should 
be noted that in general a single run of Shor's algorithm will not be sufficient because it 
is a stochastic algorithm. In the following we will discuss the time required to perform 
one run of Shor's algorithm and if not stated explicitly the calculation time is just the 
time required for this. 

The calculation time has to be compared to the decoherence time Tdec of the quantum 
computer (eg the time in which on average one photon will be emitted by the quantum 
computer). As spontaneous emissions destroy the coherence in the quantum computer, 
we need to make sure that practically no spontaneous emission occurs during the whole 
computation. To ensure this, the inequality 

Tdec > T = eTelL"^ (1) 

has to be satisfied which then gives rise to an upper limit for the numbers we are able to 
factorize on the quantum computer. For a given value of Tei that means that the total 
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computation time scales like L^. To factorize a number representable by L qubits, one 
requires 5L + 2 qubits (in what follows we neglect the "2" here) as work space for the 
necessary calculations 0. If we assume that each qubit couples to a different bath the 
decoherence time of 5L qubits is given by [|^, |[ 

r.ec = g (2) 

where r^^ the decoherence time of a single qubit. The case of qubits coupling to the same 
bath leads to smaller decoherence times T^ec Hi • Combining eq. @j and eq. @) we obtain 

Tqb > Tel 5eL^ . (3) 

Usually Tel is not assumed to be related to the decoherence time of the quantum computer. 
As we will see later this is not true in general. We will show that the dependence of 
the elementary time step t^i on the decoherence time r^ec gives rise to a much stronger 
dependence of the calculation time on the bit size L. This results in a severe limitation of 
the maximum size of the numbers to be factorized. In our investigation we focus on the 
model put forward by Cirac and ZoUer [Q but also show briefly that similar restrictions 
apply for cavity QED implementations. We stress that the results apply to a wide class 
of possible models as most of them rely on atom-light interaction similar to that of the 
models discussed here. Of course the actual form of T{L) may vary slightly from model 
to model. 

In Section II we investigate the model of a quantum computer proposed by Cirac 
and Zoller for several possible methods to store the qubits as well as a cavity QED 
implementation. In Section III we summarize our results and discuss their implications 
to the realizability of quantum computers. 

II. Quantum Computation in a linear ion trap 

In the introduction we gave a simple estimate of the time T a quantum computer requires 
to perform Shor's algorithm. From this it is possible to obtain an upper limit for the 
numbers that we are able to factorize. However in this estimate it is usually assumed that 
the execution time for an elementary logical gate does not depend on the decoherence 
time of the quantum bits on which the operations are performed. This however is not 
generally true. To see this note that all the proposals for the practical implementation of 
quantum computers mentioned in the introduction share a common feature. They rely on 
the interaction of light with atoms where either the atoms are used as a memory to store 
the qubits which are manipulated by light fields or the light field is used as the memory 
which is manipulated by the interaction with atoms. Therefore in all these schemes the 
atom-light interaction represents the essential building block of all the proposals made so 
far. In each of these interactions a temporary excitation of the atoms is inevitable (even 
in adiabatic excitation, given a finite excitation time) which can lead to spontaneous 
decay. Obviously the interaction strength, proportional to the Rabi frequency f2, and the 
spontaneous emission rate, proportional to the Einstein coefficient of the excited level of 
the transition in question, are related such that 

where F is half the Einstein coefficient of the transition and p is a constant of propor- 
tionality. Certainly for a given transition frequency p cannot be made arbitrarily large. 
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It is limited due to the fact that at high intensities the two level approximation breaks 
down, that the rotating wave approximation becomes invalid and that for a sufficiently 
high laser intensity the atom ionizes practically immediately. For optical transitions the 
latter effect gives rise to an upper limit of the order of 

Pma. = 10^°S^/' . (5) 

In practize the limit will be much lower as both detuning and pulse duration have to be 
controllable quantities and we have not included the other limitations mentioned above 
in eq. (|^) . As the execution time t^i of a quantum gate depends inversely on the Rabi fre- 
quency VL while the decoherence time of a qubit r^b depends inversely on F we immediately 
observe via eq. @) that both quantities are related to each other. 

In the following we will investigate how this relationship affects the estimate for the 
factorization time of a number which can be represented by L qubits. First we discuss the 
scheme proposed by Cirac and Zoller because it seems to be the most promising proposal. 
Later we show that for cavity QED implementations similar problems arise. In similar 
ways one may achieve estimates for other proposed schemes as they mostly rely on atom- 
light interaction. The exact form of T{L) might be different but one will always find that 
the scaling with L is much stronger than expected from eq. (|l]). 

A. Linear trap with two level atoms as qubits 

We now discuss the model proposed by Cirac and Zoller ^j. Several ions of mass M 
are stored in a linear trap (see Fig. 1) and it is assumed that all translational degrees 
of freedom of the ions are cooled to their respective ground state and that especially 
the center-of-mass (COM) motion with frequency v is in its ground state. This implies 
that the Lamb-Dicke regime is reached. To implement quantum gates one then applies a 
sequence of laser pulses of wavelength A to the ions such that both the internal degrees 
of freedom as well as the degree of excitation of the COM mode may be changed. As 
the COM mode is a collective motion of all ions, its excitation can be used to yield 
entanglement between different ions. As an approximation it is assumed that only the 
COM mode is excited because the closest lying mode has a frequency v^i/ and is therefore 
well separated from the COM mode frequency. In the model it is assumed that the laser 
is detuned such that A = —v, so that the predominant contribution comes from processes 
where with the excitation of the ion the COM mode is deexcited. Processes where the ion 
and the mode are excited simultaneously include rapidly oscillating phasefactors and are 
neglected in the following (rotating wave approximation). One then obtains the following 
Hamilton operator for an ion at the node of a standing light field 
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e)((7|a+|(7)(e|at . (6) 



where rj = ^-y/(^/2Mz/) <^ 1 is the Lamb-Dicke parameter. The a and a' are the 
annihilation and creation operators of the COM mode. The Hamiltonian eq. (|^) is 
correct for (i7/2z/)^?7^ <^ 1. This system allows the implementation of elementary logical 
gates such as the controlled-NOT gate [|l| which requires in this scheme the equivalent of 
four TT-pulses with the Hamiltonian eq. (P). We use the time required for this as a lower 
bound for the elementary time step Tei and find 
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Now using the fact that Shor's algorithm requires eL^ elementary steps we find for the 
total computation time 

4 7ry5l 



T 



7]Vt 



-eL' 



As we want to minimize T, we insert the maximum value for Q according to eq. 
obtain 
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In this expression not all the parameters are independent, as we have to make sure that 
T is less than the decoherence time r^ec of the quantum computer. The decoherence time 
of the quantum computer is the decoherence time of a single quantum bit Tgi, divided by 
the number of quantum bits contained in the quantum computer because in the course of 
the calculation most of the qubits will be partially excited. We find 
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and obtain the inequality 



We observe that due to eq. ( 
of the equation and we find 
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the decay constant of a single qubit appears on both sides 
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which is far more restrictive than the estimate eq. (|^) obtained when we assume that an 
elementary time step r^i is independent of r^gc- To be able to perform Shor's algorithm 
without having spontaneous emissions eq. (|T2p has to be satisfied. Using this to eliminate 
r in eq. (|^) then gives a lower bound for the calculation time which is 
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To obtain explicit values for T we assume rj 
the order of 1000 [01 so that we obtain 



0.1 and p = lO'^s The value of e is of 
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One observes that even with the rather large value of p the factorization of a 4 bit 
number (eg. 15 which is the smallest composite number for which Shor's algorithm 
applies [0) seems to be practically impossible when we take into account that for example 
the metastable transition in Barium has a lifetime of 45s and therefore F = 0.044s~^. 
Note that we have not taken into account the infiuence of all other possible sources of 
error such as counterrotating terms in the Hamilton operator, excitations of modes other 
than the COM mode, errors in the pulse lengths and in the Rabi frequencies of the 
pulses. One should also realize that although a heroic experimental effort might make the 
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factorization of a 4 bit number possible, the factorization of any number of relevant size 
seems completely out of question as the execution time of Shor's algorithm for a 40 bit 
number is 10^ times larger. For a 400 bit number, which represents the upper limit which 
classical computers can factorize, Shor's algorithm requires 10^^ times longer than for a 
4 bit number. 

The main problem in the model seems to be that a metastable transition cannot be 
driven very strongly which in turn severely limits the execution time of an elementary 
gate. As a possible way to improve the above model, it was proposed to consider a 
j = 1/2 j = 1/2 transition where the qubit is represented by the two lower levels of 
the transition ||^. However in the following we will show that this scheme suffers from 
similar drawbacks as the previously investigated system. 

B. The j=l/2^ j=l/2 transition 

The level scheme we now investigate is depicted in Fig. 2. A qubit is represented by 
the levels 1 and 2 which are assumed to be stable. The transition to the two upper 
levels, however, may be strong to allow for rapid transitions. As the implementation 
of quantum gates requires the excitation of one phonon in the COM mode, we need 
to transfer population between the two lower levels with a simultaneous excitation (or 
deexcitation) of the COM mode. To be able to perform this population transfer without 
appreciable population of the upper levels which would lead to spontaneous emissions. 



one has to use the method of adiabatic population transfer [|I^ . The energy levels shown 
in Fig. 3 are the most relevant. The vertical axis gives the energy of the bare states \i;n) 
where i is an atomic level and n is the number of phonons in the COM mode. Assume that 
initially the population is in level |2; 0) and we want to transfer it to level |1; 1). During 
the (quasi)-adiabatic population transfer one first applies a cr-polarized laser pulse with 
a detuning A = —u; we assume that the ion rests at the node of the light field. The 
duration of this pulse is a fixed fraction of the total length Tad of the process while the 
length Tad of the process may be varied. Later but still overlapping with the cx-polarized 
laser pulse, a pulse of vr-polarized light is applied to the same ion and it is assumed that 
the ion is situated at the antinode of this field. This pulse, in leading order, preserves 
the excitation number of the COM mode. Again its length is a certain fraction of the 
total time Tad and we assume that the cr-polarized laser pulse terminates earlier than 
the TT-polarized pulse. If the time Tad in which this process is performed is sufficiently 
long then the population in the upper level |3; 0) will be small and therefore spontaneous 
emissions rare. This method certainly has the advantage that the exact pulse shape of 
the laser is not as important as in the previously discussed scheme. At first glance it also 
appears to be possible that the population transfer can be made extremely fast as the 
Rabi frequency is not related to the lifetime of the lower levels. However there is a limit 
to the Rabi frequency. To see this we have to realize that an adiabatic process requires 
infinite time. However if we want to be able to perform the factorization in finite time we 
have to take into account small deviations from the adiabatic behaviour. In this case some 
population will end up in the excited levels which may subsequently lead to spontaneous 
emissions. We find for the probability pem that at least one spontaneous emission takes 
place during the (quasi)-adiabatic process 

Kf 1 

Pem = /^r— -— (14) 

ri^ni Tad 
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where the constant (3 depends on the peak value of the Rabi frequency Vt^^ of the vr- 
polarized laser, the pulse shapes and the delay between the pulses. Vt„ is the peak value 
of the Rabi frequency of the cr-polarized laser. If fi^ is larger than rjVta^ and F (which 
we implicitly assume in eq. (p!^) we find for sin^-pulse shapes [3 ~ 100. Analytically as 
well as numerically one finds that (3 exhibits a very slow increase with increasing VL^^. We 
have assumed that the (quasi)-adiabatic process is sufficiently slow so that the 1/T law 



applies. This is the case when the right hand side of eq. (14) is small compared to one. 
As we do not want to find any spontaneous emission during the whole computation the 
inequahty 

^-^-^=PemeL' 15 

77^ Tad 

needs to be satisfied. This gives an estimate for the length of an elementary time step 
which is 

Tel ~ Tad > ^^^^L^ ■ (16) 

Therefore we obtain for the total calculation time the estimate 
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Again this estimate scales much stronger with the bitsize L of the input than expected. 
To see the orders of magnitude, we give explicit values for T. Assuming 77 = 0.1,/? = 
100, e = 1000 and p = Ws'^^^ we obtain 
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which indicates that even the factorization of a 4 bit number will be extremely difficult to 
achieve, although the estimate seems to be a little more promising than in the previous 
scheme. Again we have neglected all other sources of error, such as higher order contribu- 
tions in the Lamb-Dicke parameter to the Hamilton operator as well as counterrotating 
contributions neglected in the rotating wave approximation. Because the expression eq. 
(|17|) contains the ratio T/Q"^, again we have similar problems as before as this ratio cannot 
be made arbitrarily small. 

C. Cavity QED implementation 

Now we would like to show briefiy that in cavity QED realizations of quantum computing 
expressions similar to eq. (|13D and eq. ([T7|) can be obtained. In cavity QED implemen- 
tations of quantum gates the atom-light interaction does not involve a classical laser field 
but a quantized mode of a cavity. Before and after the cavity we may use Ramsey zones 
to rotate the Bloch vector of the atoms passing the cavity 0|. To perform quantum com- 
putations such as Shor's algorithm, many cavities are required and this obviously poses 
immense experimental difficulties. In the following we neglect the restrictions arising from 
these problems as well as all difficulties that arise in the realization of exactly one atom 
passing with a well defined velocity through the cavity. We will briefiy show that again 
the lower bound for the computation time scales much stronger than with the bit size 
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L of the number to be factorized. Neglecting decay of the cavity mode, we can estimate 
that the minimal computation time is of the order of 



Train ^ (-^^) 

where Vt is the Rabi frequency in the cavity-atom interaction. While travelling in the 
Ramsey zones and between cavities the atoms may decay. No decay should occur during 
the quantum computationcal which leads to the condition 

"-La^ « 1 («) 

where a depends on the ratio between the time the ion spends inside the cavity (where we 
neglect spontaneous decay) to the time it spends outside the cavity (where it may decay). 
Using eq. (^) we then obtain 

T > . 20 

Although this estimate seems much more promising than eq. (p!3[) and eq. (|l3), it should 
be noted that it is certainly an unrealistically low limit because we have negelcted major 
sources of experimental uncertainty mentioned above. We only intend to illustrate that 
again an expression similar to eq. (|TB|) and eq. (0) is found although we have discussed 
a completely different realization. 



These examples show that it seems to be a general feature that the control of pop- 
ulation always leads to the appearance of a factor of the form F/fi^ which, for a given 
transition frequency, has an upper limit. There seems to be only one way out of this 
dilemma. Instead of employing optical transitions to represent qubits one could use low 
frequency transitions (e.g. microwave transitions) as it was done in the cavity QED imple- 
mentation of Sleator and Weinfurter 0] because this can considerably decrease the ratio 
r/Jl^ = Xj (? due to the ur" dependence of F. However as in their proposal one would 
need a tremendous number of cavities it does not seem very promising. To overcome 
this problem one might use the cavity field in the manner implementation by Cirac and 
ZoUer 0. Instead of using the COM mode to entangle different ions this task could be 
performed by the cavity mode. This could be done using a linear trap to store the ions 
inside a microwave cavity. This scheme then resembles that of Sleator and Weinfurter 
but differs as we only require one cavity and we do not need atomic beams with all their 
associated problems. The COM mode will not be excited during the calculation as for the 
long wavelength of the radiation the Lamb-Dicke parameter is extremely small. However 
smaller frequencies of the incident fields mean larger wavelengths which will make it more 
difficult to address single ions with the microwave radiation. The problem of addressing 
a single ion, given many are within a wavelength of the incident radiation, may be solved 
by applying local magnetic or electric fields (or a suitable field gradient) that drive all 
but one ion out of resonance. However due to the small spatial separation of the ions this 
might be difficult to realize experimentally. If it would be possible to implement this idea 
then the lowest limit for the computation time could become as low as eq. (^) with a 
value of p that can be much larger than that for an optical transition. However this idea 
should serve rather as a basis for discussions than a serious proposal as we still expect 
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the experimental difficulties to be enormous. We are therefore not very optimistic that 
factorization of nontrivial numbers will be possible in the near future. 

III. Summary 

In this paper we have investigated how the computation time which a quantum computer 
needs to factorize an L bit number depends on several physical parameters. It was shown 
that T will scale much stronger with L than previously expected. Instead of an Li^ de- 
pendence we find an or behaviour in the proposal of Cirac and ZoUer and for 
cavity QED realizations in which however this limit is more of theoretical nature than of 
practical importance due to other experimental problems. In the models that we have 
investigated explicitly, it also turns out that the computation time is always dependent 
on the ratio T /VL^ where F and are the decay constant and the Rabi frequency of one 
of the transitions that are required to transfer population. Although found for special 
configurations, this seems to be a general result which limits the length of the elementary 
time step because the ratio T /Vt^ cannot be made arbitrarily small for an optical tran- 
sition. As a possible way to circumvent these problems, we briefly discussed the use of 
microwave transitions to store qubits as in this case the ratio V/Vl? becomes extremely 
small. However practical problems occur which seem to make the experimental realiza- 
tion of this idea difficult, although it might lead at least to the possibility to factorize 
numbers which are several bits long, a task which seems to be impossible with the present 
proposals. 
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FIGURE CAPTIONS 



Fig. 1 : Schematic picture of the excitation of several ions in a hnear ion trap. The trans- 
lational degrees of freedom of the ions are assumed to be cooled to their respective 
ground states. To implement quantum gates, standing wave fields interact with 
the ions and thereby changing the inner state of the ions as well as the state of 
the center-of-mass mode which leads to entanglement. 

Fig. 2: A J = 1/2^ J = 1/2 transition. The qubit is represented by the two lower levels 
1 and 2. Population transfer requires two different lasers. Adiabatic population 
transfer minimizes unwanted population in the upper level. 

Fig. 3 : The j = 1/2 ^ j = 1/2 transition including the quantized center-of-mass mo- 
tion. \i; n) denotes an atomic level i and n phonon in the center-of-mass mode. 
For the implementation of a controUed-NOT gate we need to be able to transfer 
population from state |2; 0) to state |1; 1) and vice versa. To minimize population 
in the excited levels population transfer is performed using adiabatic population 
transfer with counterintuitive pulse sequence. 
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Standing wave laser fields 
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